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We report measurements of thermally- induced mean-square director-angle fluctuations (8 2 ) below 
primary supercritical bifurcations to electroconvection of the nematic liquid crystals 152 and "Merck 
Phase V". For e m f = V 2 /V 2 tUlf — 1 < —0.1 (V is the amplitude of the applied alternating voltage) 
we find (8 2 ) oc |e m f| ' with 7 given by linear theory (LT). Closer to the bifurcation there are 
deviations from LT with a smaller 7 and with V 2 > V^. 2 mf . For 152 measurements as a function of 
the conductivity a of {8 2 ) above V 2 suggest a tricritical bifurcation at at — 4.0 x W~ 9 Q,~ 1 m~ 1 , and 
(V 2 — K. 2 mf )/V c 2 increases strongly as a t is approached from above. 
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Convection in fluids is a classical system for the study 
of pattern formation under noncquilibrium conditions 
because of its similarity to numerous pattern-forming 
phenomena in nature Q. Rayleigh-Benard convection 
(RBC) is the most prominent example. Here a horizontal 
fluid layer, confined at the bottom and top, is heated from 
below. If the temperature difference exceeds a threshold, 
the system undergoes a sharp but continuous (supercriti- 
cal) bifurcation from a uniform state to a state where pat- 
terns, e.g., stripes, squares, or spirals, occur ||. However, 
already below this bifurcation there exist fluctuations ST 
of the temperature field which become "large" near the 
bifurcation point and which are induced by thermal noise. 
The fluctuation amplitudes have zero mean but a finite 
mean-square (ST 2 ). For RBC quantitative predictions of 
(ST 2 ) were made three decades ago [§|,|],|]|6| on the basis 
of linear theory (LT) which neglects interactions between 
the fluctuations. Since (ST 2 ) is extremely small, quanti- 
tative experimental verifications could be obtained only 
much more recently . 

Closer to the bifurcation, LT should break down be- 
cause of nonlinear interactions between the fluctuations. 
In analogy to critical phenomena in equilibrium systems, 
one then expects a modified "critical" (rather than "mean 
field") behavior of the system. For RBC it was predicted 
by Swift and Hohenberg || that the fluctuation inter- 
actions should lead to a first-order transition, i.e., to a 
subcritical bifurcation. Experimentally this interesting 
phenomenon has been out of reach so far because under 
most circumstances it is expected to become noticeable 
only within a few parts per million of the bifurcation 
point 



with NLC. The role of the temperature difference is now 
taken by V, and EC occurs for V > V c . For NLCs the ef- 
fects of thermal noise are larger than for RBC because of 
the small elastic constants of the NLC pT| . This already 
large susceptibility is enhanced even further in EC by the 
very small thickness of the cell. Indeed, already a decade 
ago it was possible to visualize the fluctuating convective 
patches and to measure their amplitudes below onset [fl2| . 
Good agreement with LT was found. We do not know of 
a detailed nonlinear theory for this system. Its critical 
behavior need not be the same as that of RBC because 
the anisotropy of the NLC suggests that it may belong to 
a different universality class. Although one expects the 
critical region to be wider than for RBC, its study had 
been beyond experimental resolution heretofore. 

We report results for thermally-induced fluctuations of 
the director angle 9 of the NLCs 152 and "Merck Phase 
V" (MPV) both far below and very close to the primary 
supercritical bifurcation. For 152 we measured the mean 
square (0 2 ), averaged over time as well as over all fluctu- 



A more favorable system to study the influence of ther- 
mal noise in pattern-forming systems is electroconvection 
(EC) in a nematic liquid crystal (NLC) p0[ . Here an al- 
ternating voltage of amplitude V is applied to a cell filled 



ating modes of the system as described in [ |13|jl4[| . In the 
case of MPV, however, we measured the time averaged 
mean square amplitudes (A 2 ) of the fluctuations of only 
the x-component of the critical mode of the system, as de- 
scribed in |ll| . In agreement with earlier work, |l^,^3|,Q 
our data well below V c are consistent with (6 2 ) oc e~f 
((A 2 ) cx e~J°) with 7 = 1/2 (70 = 1) as predicted by 
LT for the multi-mode (single-mode) measurements (here 
e m f = V 2 /V 2 mf — 1). The data extrapolate to a mean- 
field threshold at T4 jm f. However, closer to the bifurca- 
tion we find (6 2 ) oc eT 1 with e = V 2 /V 2 — 1, a smaller 
7 ~ 0.22 (70 ~ 0.5) for the multi-mode (single-mode) 
measurements, and a shifted threshold V c > V^mf. So 
far as we can tell the bifurcation remains supercritical 
unlike the prediction for RBC. We suggest that an ex- 
planation of this deviation from LT may be found in the 
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nonlinear interactions of the fluctuations. Unfortunately 
there is as yet no detailed theory for the critical behavior 
of this system with which our results could be compared. 

TABLE I. Parameters and results for 152 at different 
temperatures i9. 
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FIG. 1. Background-divided and Fourier-filtered snap- 
shots of size 502 x 502 (im 2 (a,c), and the central part 
— 22.4/d < k x , k y < 22.4/d of the time-averaged structure fac- 
tors 5(k) (b,d) for 152. (a,b) : e = -0.104. (c,d) : e = -0.001. 
The director is horizontal. 

In the experiments reported here we used samples with 
planar alignment. The 152 (4-ethyl-2-fluoro-4' -[2-(trans- 
4-pentylcyclohexyl)ethyl]-biphenyl) was doped with 4.3% 
by weight of molecular iodine (±2) and studied in the ap- 
paratus and by the cell-assembly techniques of Ref . [Q . 
The bifurcation to EC is supercritical over a range of the 
electrical conductivity a . Above the bifurcation the 
pattern consists of chaotic traveling zig and zag modes. 
The cell had a thickness d = 28 ± 1/xm. We changed a 
by varying the temperature (see Tab. and measured 
a at a frequency of 50 Hz and F = 2.0 V. For the ex- 
periment we used a frequency of 25 Hz. At each voltage 
we waited 110 s, and then took 128 images 10 s apart. 
The combination of cell thickness, I2 concentration, and 
drive frequency assured that the "worm" state (l7| did 
not occur for a > 4 x 10 
frequency of 30 Hz and d 

Figure [l] shows two examples of single snapshots for 



9 fi _1 m -:i . For MPV we used a 
23.1 ±0.2 ^m HI. 



152 at different e and for a = 7.8 x 10 _9 O _1 m _1 . Far 
below onset [Fig. 0(a)] we see patches of patterns with 
a small correlation length. The patterns are very weak, 
although dividing by a background image, Fourier filter- 
ing, and using the full available greyscale uncovers the 
expected zig and zag modes. Much closer to onset [Fig. 
|l|(c)] a pattern containing extended patches of zig and 
zag rolls is found. 
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FIG. 2. (a): The root-mean-square of the director-angle 
fluctuations (O 2 ) 1 ^ 2 and (b): the inverse square of the mean 
square of the director-angle fluctuations (8 2 }~ 2 for 152 and 
a = 7.8 x 10 _9 n _1 m _1 as a function of the square of the 
driving voltage. The dotted lines are a fit of a straight line 
to (6 2 )(V 2 ) in the range 124.53 V 2 < V 2 < 125.15 V 2 . 
They extrapolate to V 2 = 124.495 V 2 . The dashed curves 
are a fit of a straight line to the data in (b) in the region 
85 V 2 < V 2 < 114 V 2 where LT applies. In (b) the inter- 
cept of the dashed straight line with the horizontal axis gives 



For the multi-mode analysis with 152 we closely fol- 



lowed the method of Ref. 



ilysi; 




to extract (9 ). We cal- 



culated 7i(x, e) = Ii (x, e)/7o (x, e) — 1 for each image 
Ji(x,e), i = 1,...,128. Here x = (x,y) are the coordi- 
nates in real space and Io(x,e) is a background image 
obtained by averaging 128 images at the same e. For 
each Ii (x, e) we derived the structure factor (the square 
of the modulus of the Fourier transform) Si (k, e) and av- 
eraged 128 Si(k,e) to get 5(k, e), where k = (k x ,k y ) is 
the wave vector. Figures JL|b) and (d) are examples. In 
agreement with Refs. [|l3],[t4 jl8 1 we see two pairs of peaks 
corresponding to two sets of rolls oriented obliquely to 
the director which get sharper and larger as we approach 
the onset of convection. The two modes are called zig 
and zag modes and correspond to those of the extended 
chaos above onset |lq|. We computed the total power 
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under the peaks of 5(k) as described in [pTl, and then 
converted it to (8 2 ) (ij (see also Eq. (5) inflj]). 

Figure ||(a) shows (9 2 ) 1 / 2 as a function of V 2 for 152 
and a = 7.8 x 10~ 9 f2 m -1 . It is small for low voltages 
but then increases sharply above onset. As expected, 
the primary bifurcation is supercritical in this interme- 
diate region of a, and the points above onset follow a 
square-root law. The square-root fit above onset [dotted 
line in (a)] extrapolates to (0 2 ) 1 / 2 = at V 2 . 
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FIG. 3. (a): Mean square (9 2 ) for 152, averaged over all 
modes and over 128 images, of the director angle as a func- 
tion of log ( | e | ) (open circles: e < 0, solid squares: e > 0). 
Short dashed curve: fit to data in the linear region from 
Fig. |^(b). Straight solid line: fit for e < in the range 
-3.31 < log(|e|) < -2.13. It yields 7 = 0.22. Dotted line: 
fit to data above onset from Fig. ^(a). Vertical dashed line: 
location of V c , m { along the e-axis. (b): Mean square {A 2 } of 
single-mode fluctuation amplitudes (arbitrary scale) for MPV 
as a function of log(|e|). The lines and symbols have meanings 
equivalent to those in (a), (c): Effective local exponent 7 e ff 
of (e 2 ) in (a) for e < 0. It was obtained from fits within a 
sliding window covering half a decade of e. Dashed horizontal 
line: 7 = 0.5. Solid horizontal line: 7 = 0.22 as suggested by 
our data at small |e|. 

In Fig. |(b) we plotted (6> 2 }~ 2 versus V 2 . The data 
are the same as those in Fig. |^(a). For V 2 < 113.6 V 2 
they can be fitted by a straight line. This is in agree- 
ment with LT, which gives (6 2 ) = 9 2 ,e^ 2 . The straight- 
line fit shown in the linear region of Fig. |(b) yielded 
V 2 m{ = 122.685 V 2 and 9 = 1.8 mrad. The value of 
9q is in fairly good agreement with the approximate the- 
oretical estimate 9q = 2.8 mrad (see Ref. |l2| and Eq. 



(6) in |14|). Hence, the measured mean-square director- 
angle fluctuations are consistent with thermal noise. The 
mean-field threshold V 2 mf is well below the actual bifur- 
cation point at V^ 2 . A similar, albeit somewhat smaller, 
shift was found for MPV (see Fig. |(b) below). 

Figure |](b) reveals that the points close to but below 
onset deviate from the prediction of LT. Similar devia- 
tions were observed over a range of a (see Tab. |), as well 
as for MPV. To illustrate the behavior close to the bi- 
furcation, we plotted log((0 2 )) versus log(|e|) in Fig. ft(a) 
for 152 and a = 7.8 x 10 -9 fi -1 m -1 . For e < (open 
circles) the data points far away from V c (far right of 
the figure) follow the behavior predicted by LT (short- 
dashed curve). Note that in a plot where the horizontal 
axis is log(|e m f|) this should give a straight line with a 
slope of -1/2 for the multi-mode analysis used here. The 
line is curved and diverges at the e-value corresponding 
to T^, 2 mf . Closer to onset the data points deviate from 
the linear prediction. We find a crossover to a regime 
which can be described by a new powerlaw: (9 2 ) oc e -7 
with 7 ~ 0.22. In order to characterize this crossover, 
we calculated the local exponent 7 c ff (e) from fits to data 
spanning half a decade of e. Figure 0(c) shows that 7 e g 
matches the expected 7 = 0.5 in the linear region at large 
|e|. As we approach the onset of convection, the exponent 
approaches a plateau of 0.22. Thus we plotted a line in 
Fig. ||(a) with the slope of -0.22 (straight solid line) and 
found good agreement with the data over a wide range. 
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FIG. 4. Mean square (6 ) of the director angle 8 close 
to onset as a function of e for different a. Straight dotted 
lines are fits above onset. Open circles, solid circles, open 
squares, and solid squares correspond to a — 5.2, 6.4, 7.8, and 
9.5 x 10 _9 fl _1 m _1 , respectively. 

Data for MPV are shown for comparison in Fig. ||(b). 
As mentioned above, they are based on a single-mode 
analysis. Far below the bifurcation they yield the ex- 
pected LT exponent 70 = 1 (note that the resolutions of 
the vertical scales of Figs. ||(a) and (b) differ by a factor 
of two so that the corresponding slopes look equal) . With 
decreasing |e| there is a crossover to a smaller exponent 
close to 70 = 0.5, which is consistent with the crossover 
from 7 = 0.5 to 7 = 0.22 found from the multimode- 
analysis for 152. 

For 152, measurements similar to those discussed above 
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were made at four different conductivities (see Tab. |). 
Figure || shows the results in the vicinity of onset. The 
data for e > can be fit by (8 2 ) = Si(a)e as expected 
for a supercritical bifurcation. The slope Si increases 
as a decreases. In Fig. || we show l/Si (open circles) 
as a function of a. The data suggest that \/S\ van- 
ishes at at ~ 4.0 x 10 _9 r2 _1 m _1 . This is consistent with 
a tricritical bifurcation at a t , with a transition from a 
supercritical to a subcritical bifurcation as a decreases 
below at- The results for 1 /Si are consistent with recent 
weakly- nonlinear calculations pOf . 
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FIG. 5. The inverses of the initial slope Si (open circles) of 
(9 2 ) above onset, and of the onset shift Ae (solid squares: this 
work, open square: from Ref. jL4|) as a function of a. The 
solid line is a fit to 1/Si. Its intercept with 1/Si = suggests 
a tricritical point at at = 4.0 x 10 _9 Jl _1 m _1 . A straight-line 
fit to 1/Ae also passes through zero near at- 



In Fig. H we also show 1/Ae where Ae = (V c — 
V C}ln s)/V c . Here the solid squares are from the present 
work, and the open square is from @ The dashed 
straight line is a least-squares fit to the data. Within our 
resolution it passes through zero at o~ t , suggesting that 
the critical region broadens strongly near at |22|. 



In this paper we showed that the fluctuations below 
onset of electroconvection reveal a crossover as the bifur- 
cation is approached from the behavior predicted by lin- 
ear theory to a different region characterized by a smaller 
exponent and a shifted threshold. In analogy to equilib- 
rium phase transitions, we suggest that this phenomenon 
may be attributable to nonlinear interactions between the 
fluctuations. Equivalent behavior was found in the two 
NLCs 152 and MPV, even though we used very different 
methods of analysis for the two cases. 
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